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1 Introduction 


Entanglement is a fundamental property of quantum systems which relates to the outcomes of 
local measurements. In the context of 1+1 dimensional many body quantum systems, especially 
over the last 10 years, there has been growing interest in developing efficient (theoretical) mea¬ 
sures of entanglement as a means to extract valuable information about emergent properties of 
quantum states. Various measures of entanglement exist, see e.g. [1, 2, 3, 4, 5, 6, 7], which occur 
in the context of quantum computing. Among these, a well-known and much studied measure 
of entanglement within many-body physics is the bipartite entanglement entropy [1] (see [8] 
for various reviews). Another measure that has recently received some attention within many- 
body physics [9, 10] is the logarithmic negativity [6], which was shown to be an entanglement 
monotone in [7]. 

The entanglement entropy (EE) measures the amount of quantum entanglement, in a pure 
quantum state, between the degrees of freedom associated to two sets of independent observables 
whose union is a complete set for the Hilbert space. Computations of the EE have revealed them¬ 
selves as powerful tools especially in the study of quantum critical points, since the EE exhibits 
striking universal behaviours at and near such points. In 1+1 dimensions, the behaviour at crit¬ 
icality has been established analytically by employing conformal held theory (CFT) techniques, 
which have proven extremely effective in predicting the EE scaling laws in spatial bipartitions. 
The first of such results [11, 12] established a fundamental relation between the logarithimic 
growth of the EE and the CFT central charge for connected regions in ground states of unitary 
models, later re-derived using different techniques [13, 14]. Generalizations have included ex¬ 
cited states [15, 16], disconnected regions [17, 18, 19, 20, 21] and non-unitary models [22], The 
first result beyond criticality [13] showed that the EE of infinite regions saturates to a value 
which logarithmically grows with the correlation length. In the full universal region near critical 
points, the EE is described by a universal scaling function related to correlation functions of 
branch point twist fields in massive quantum field theory (QFT) [23]. Using this, a precise de¬ 
scription of corrections to the saturation value has been found [23, 24], showing that the leading 
system-size dependent correction displays highly universal and unexpected features controlled 
by the mass spectrum of the QFT. 
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Figure 1: Typical negativity configurations in 1+1- and 1+2-dimensions. 

The EE, however, is not a good measure of entanglement in mixed states. This, for instance, 
precludes its use for measuring the entanglement between non-complementary regions, such 
as the regions A and B depicted in Fig. 1, since the reduced density matrix on the subspace 
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corresponding to A U B is generically mixed. An alternative measure of entanglement which is 
suited to treat mixed states is the logarithmic negativity introduced in [6, 7]. Powerful methods 
were developed in [9, 10] which allow its study in CFT using branch-point twist fields, much like 
in the case of the EE, and some universal results were found for critical points, again relating 
it to the CFT central charge. However, still relatively little is known about other universal 
features of the negativity. Detailed studies of free massless theories such as the compactified 
free Boson [10, 25], critical Ising spin chain [26] and free massless Fermion [27] have been carried 
out, as well as investigations of the negativity after a quantum quench [28, 29, 30, 31] and at 
finite temperature [32]. The negativity has also been recently considered within the holographic 
approach [33, 34], 

In this paper, we study the logarithmic negativity in general unitary 1+1 dimensional massive 
QFT. We consider two cases of the configuration presented in Fig. 1: (i) the regions A and B 
are adjacent, B is semi-infinite and A covers a finite but large distance; and (ii) both A and 
B are semi-infinite and separated by a finite but large distance. We obtain universal results 
for large-distance asymptotics purely in terms of QFT data. In particular, we find, like in the 
case of the EE, that the leading correction to the large distance constant (there is large-distance 
saturation in case (i), and vanishing in case (ii)) is completely determined by the mass spectrum 
of the model. Interestingly, unlike for the EE, we find that certain bound state data are involved 
in the first few sub-leading terms of the large-distance expansion in case (i). 

This paper is organized as follows: in section 2 we provide a formal definition of the negativity 
and review its representation in terms of a four-point function of branch point twist fields [9, 10]. 
In section 3 we describe in- and out- states and the analytic properties of the matrix elements 
of branch-point twist fields in general massive 1+1 dimensional (replica) QFT [24], In section 4 
we present a completely general computation of the leading correction to negativity saturation 
in massive 1+1 dimensional QFT for the cases of two adjacent and of two semi-infinite non- 
adjacent regions. In section 5 we present our conclusions and outlook. Specific examples are 
briefly discussed in appendix A. 


2 Logarithmic negativity in quantum field theory 


The logarithmic negativity provides a measure of the bipartite entanglement between two sets 
of observables in mixed states. Consider the geometric configuration presented in Fig. 1, with 
the whole system being in a pure state Tracing over the space corresponding to C, the 
logarithmic negativity of the resulting reduced density matrix provides a measure of the bipartite 
entanglement between regions A and B. The precise definition goes as follows: let p be the 
reduced density matrix p = Trc(|'i/’)(V>|) of the subsystem A L) B and consider the “partial 
transpose” of this density matrix, denoted by p Te , with respect to the tensor factor corresponding 
to B. That is, let \ef) and \e^), for i,j e N, be basis state vectors on the Hilbert spaces 
associated to the degrees of freedom on A and B respectively, and denote = \e^)A® |e ?)b. 

The matrix elements of p Ts are then given by the relations (efe^\p TB \e^ef) = {efef\p\e^e^). 
Although this notion is basis-dependent (because the notion of transpose is not fundamental to 
a Hilbert space), it was shown in [7] that the logarithmic negativity, defined as 

£ = log Trip 7 ® | (1) 
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where Tr|p Ts | represents the sum of the absolute values of the eigenvalues of p Te , is a good 
measure of entanglement in mixed states. 

As is the case for the EE, QFT techniques are particularly well suited to the study of the 
negativity by using a replica approach. Such approach was described in detail in [9, 10] for 
CFT and is easily generalizable to unitary massive QFT and quantum chains. The description 
provided for CFT applies directly to the case we want to study. The idea is to compute 

£[n] := logT'r(/) Ts ) ri (2) 

for positive integer values of n, and to “analytically continue” the resulting function of n. A key 
conclusion from the analysis presented in [9] is that the expression (2) has different “natural” 
analytic continuations £ e [n] and £ 0 [n] for different parities of n (even and odd, respectively), 
and that the negativity (1) is obtained by taking the limit n —> 1 from the analytic function 
£ e [n] in the even sector. That is: 

£ = lim £ e [n]. (3) 

n—> 1 

It was then found that Tr {p Ts ) n can be expressed, for integer values of n , in terms of correlation 
functions of local fields in a replica theory, composed of n independant copies of the original 
model. The fields involved are the branch-point twist fields [35, 36, 13, 23], and the techniques 
and expressions obtained are similar to those used in the computation of the entanglement 
entropy [13, 14, 23]. 

Branch-point twist fields are twist fields associated with elements of the permutation symme¬ 
try group of any replica (n-copy) model. Labelling the copies by j = 1 ,,n, the branch-point 
twist held T is associated with the cyclic permutation 1 j h > j + 1 modn, and the “anti-twist 
held” T is associated with the opposice permutation j i—> j — 1 mod n (these generate permu¬ 
tation cycles of length n ). The locality properties of branch-point twist helds were described in 
detail in [23], leading to a set of consistency equations for their matrix elements in integrable 
1+1-dimensional QFT, a part of which still holds in arbitrary massive QFT [24]. In CFT, the 
conformal dimensions of both twist helds are known to be [35, 36, 13] 



where c is the central charge of the CFT. In general, this formula holds for twist helds, in 
replica models, associated to any permutation element, where n is the length of the permutation 
cycle generated. We normalize branch-point twist helds in the usual conformal normalization, 
represented by the OPE 

T(0)T(r) ~ r~ 4An l + .... (5) 

Under this normalization, the CFT structure constant Cj-tj- is unambiguous, 

T(0)T(r)~C^T 2 (0) + ... (6) 

where T 2 is the composite twist held associated with the permutation element j j + 2 mod n. 
The structure constant C-yj- will play a role below. 

In terms of these helds, the negativity associated to a configuration such as that represented 
in Fig. 1 may be obtained in terms of a four point function of twist helds located at the boundary 

1 Our convention is that the twist field T has a branch cut to its right, and the continuity condition along this 

branch cut connects sheet number j + 1 below the cut to sheet number j above it. 
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points of regions A and B [9]. Denoting the regions A = [r 4 , 7 - 2 ] and B = [ 7 * 3 , 7 - 4 ] with r± < r 2 < 
r 3 < r 4 , the quantity £[n] can be written in terms of the following four-points correlation function 

£[n\ = 8 A„loge + log (n(vac|T(ri)T(r 2 )T(r 3 )T(r 4 )|vac) n ^ (7) 

where |vac) n is the n-tensor product of the vaccum state of the original theory, |vac) n = |vac) (g> 
• • • (g) |vac). Here e is a non-universal short-distance cutoff, and the second term is a universal 
scaling function. 

The correlation function above can be very involved, even for CFTs (see [10] for a detailed 
study of the compactified massless free Boson), so that obtaining the analytic continuation £ e [n] 
and taking the limit n —> 1 remains challenging but for the simplest of cases. In massive QFTs, 
even integrable ones, the evaluation of the four point function and its analytic continuation will 
be challenging. However, there are two special limiting cases which can be studied in much 
generality. We will consider those in this paper. 


2.1 First case: adjacent regions 

The first case, corresponds to the negativity of adjacent regions , one of the regions being 
semi-infinite : the configuration for which r 3 — > r 2 '■= r and r 4 —> 00 (we will choose r\ =0). We 
will show that, in this case, the negativity behaves as follows: 

C _L mir—^O c . . . 

£ ~ -log(r/e) 

m = >1 log(mie) + £ sat - —^L- ^2 K 0 (V3m a r) + 0 {e~ Zmir ) (8) 

^ OL 

where Z > y/2> is some number that depends on the presence of particle-anti-particle bound 
states in the QFT model (if there are no such bound states, then Z = 2). Here e is a non- 
universal short-distance cutoff which may be conveniently chosen so as to ensure there are no 
0(1) corrections to the first line of ( 8 ), {m a } is the mass spectrum of the QFT with mi the 
lightest mass, and Ko(x) is the modified Bessel function. The universal saturation constant is 
then given by 

£ sat = 2 log (mf 1 (vac|T|vac)i) - log(Ci), (9) 

where Ci is the CFT structure constant Cj-fi analytically continued from even number of copies 
n to n = 1 , and where the expectation value is likewise understood as an analytic continuation 
in the number of copies. We remark that the saturation constant involves both a massive QFT 
quantity, the vacuum expectation value, and a pure CFT quantity, the structure constant Cj. 

The limit r 3 —> r 2 was studied in CFT in [9, 10] and is fully understood in this context. 
This limit is singular on the universal scaling function (as the fields T{r 2 ) and T(r 3 ) collide), a 
signal that the leading non-universal divergency is modified. The limit is of course well defined 
on the negativity itself (which has a non-universal part controlled by the miscrocopic theory), 
and after taking this limit, there is a new universal scaling function. This new scaling function 
can be obtained by using the composite field T 2 (r), associated to the permutation element 
j 1 '? j + 2 modn, in place of the product of fields T(r’ 2 )T(r 3 ). It is simple to see that for 
even values of n, this permutation element factorizes into permutation elements on even- and 
odd-numbered copies, each generating cycles of length nj2 . As a consequence, the composite 
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twist fields factorizes, for even n, as a product of twist fields acting on even-numbered and 
odd-numbered sheets: 

T 2 = T\ j even ® Vj odd ( n even ) ( 10 ) 

In particular, for n even its conformal dimension is 

A-j- 2 (n) = 2A| (n even). (11) 

Note that by opposition, T 2 corresponds to a (non-sequential) cyclic permutation generating 
a full cycle of length n when n is odd. The field T 2 is, as usual, associated to the opposite 
permutation. In (7), the dimension of the field T 2 for n even implies that, in the configuration 
where r 2 = ry, the non-universal divergence of 8 e [n] is instead of the form 4 ^A n + A»j loge, 
where e is a (possibly different) short-distance cutoff. 

We will further take the limit r 4 —> 00 , so that one of the regions is of infinite length. In 
massive QFT, correlation functions factorize: 

lim n(vac|T(ri)T 2 (r)T(r 4 )|vac) n = n (vac|T(ri)T 2 (r)|vac) n n (vac|T|vac) n . 

mir^—>00 

This means that the negativity saturates, as r 4 —> 00 , to a value that depends on r — r 4 : the 
quantity of entanglement between a semi-infinite region and a finite region is finite thanks to 
the finite correlation length. Here, the additional term log ( n (vac|T(r 4 )|vac) n ) can be absorbed 
into a re-definition of the short-distance cutoff e. The negativity is then related to the function 
(taking r\ = 0) 

£e[n] = 4 ^A n + A|j loge + log (^ n (vac|7'(0)T 2 (r)|vac) n ^ (12) 

where the expectation value is analytically continued from positive even values of n. The uni¬ 
versal information is encoded into the two-point function, and an analysis of the above, using 
Ai = — j|t, will show that this leads to (8) (upon a re-definition of e). 

2.2 Second case: disjoint semi-infinite intervals 

Another interesting limiting case corresponds to the negativity of two semi-infinite non- 
adjacent regions, namely r 4 —> — 00 , r 4 —> 00 , and we will choose r 2 = 0, r 3 = r. In this case, by 
factorization of correlation functions, the function h [n] takes the form 

£e h M = 8A n loge + log ( n (vac|T(0)T(r)|vac) n j , (13) 

where e is some non-universal short-distance cutoff and again the correlation function is analyt¬ 
ically continued from positive even values of n. We will show that this leads to the following: 

£ Hh mir-+o _|i og ( mir .) +£ shift 

1 = >>1 ^2 m “ r [ m arKo(m a r ) 2 + A' 0 (m a r)/ii(m.ar) - m a rifi(m a r) 2 ] + 0 (e ~ 4mir ) 

a 

(14) 

where the shift constant is 

£shift = 2 log (mf 1 (vac|T|vac) 1 ) + log(C'i). (15) 
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As before, C\ is the CFT structure constant Cj^j- analytically continued to n = 1 from even copy 
numbers, {m a } is the mass spectrum of the QFT with mi the smallest mass, and A'o(x), K\(x) 
are modified Bessel functions. Note that both short- and large-distance behaviours are universal, 
and that the shift constant £ s hift is closely related to the saturation constant £ sa t given in (9). 
In particular, it is possible to extract from these universal quantities both the massive QFT and 
the pure CFT parts: 

2 log (mf ^ (vac|T|vac) ■) = £,hltt + £a “ , log(Ci) = £ ’ ,h * ~ £ui . (16) 


3 Scattering theory and form factors 


Let us recall elements of the general scattering theory of massive QFT which will allow us to 
extract universal results for the negativity. 

The space of states of massive QFT has two natural bases: the in-states and the out- 
states. These are scattering states, representing asymptotically free particles either in the far 
past or the far future. Both are labelled by the rapidity of the particles, and we denote by 
\9\ ■ ■ ■ the state with asymptotic particles of rapidities ..., 9e and internal quantum 

numbers aq,... , aq, with the convention that it is an “in” state if 9\ > • • • > #£, and an “out” 
states in the opposite order. The conjugate is denoted as \9\ - • • 9()' a i r . w = (9i---0 e \. 

The analytic continuation, in the space of rapidities, of (normalized) matrix elements of 
local fields of the form (vac| 0 | 0 i • • • / (vac| 0 |vac) are referred to as form factors and 

denoted F® i a (9 1 ■ ■ ■ 9(). These are in general complicated functions of the rapidities, but they 
possess certain analytic properties [37, 38, 39, 40]. The two-particle form factors (<?i, # 2 ) 

are particularly simple: by relativistic invariance, for spinless fields they only depend on the 
rapidity difference 9 = 9\ — 62 , and as functions of this variable, they are analytic in the physical 
strip lm(0) E [0, 2-k] except for possible bound state poles on Re(0) = 0, lm(0) E (0, 27r); the 
boundary values on ( 0 , 00) and on (27n' — 00, 2iri) are in-state matrix elements, and those on 
(—oo,0) and (27ri, 27ri + 00) are out-state matrix elements [37]. 

Branch-point twist fields were studied within the general scattering theory in [24] , following 
their study [23] in integrable QFT. Recall that this twist field is defined on an n-copy replica 
model, whose Hilbert space is the n-times tensor power 'H <2n of the original Hilbert space H. of 
the QFT. In such a replica model, any internal quantum number /u = (ct, j) is labeled both by 
a quantum number a of the underlying QFT model and a copy number j £ {1,..., n}. States 
will be denoted by \ 0 \ ■ • • and since the n copies are independent, states are direct 

products, for instance 

l^l 6 l 2 )(a 1 ,l) > (Q 2 , 2 );n = l#l)ai ® |# 2 )a 2 ® |vac) < 8 > • • - <S> [vac) (17) 

where there are n — 2 vacuum factors. 


Branch-point twist fields, however, are not factorizable in this way and therefore their matrix 
elements are generally non-trivial. It was found [23, 24] that the two-particle form factors of the 
branch-point twist field have an extended physical strip. That is, the function 


F r 

Ati,/22;n 


(9 1 - 02) : = 


n { Vac | T| 9 1 02 ) Hi \n 
„(vac|T|vac) n 


(18) 
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defined as the analytic continuation of the matrix elements in terms of the variable 9 = 9 1 — 02, 
is analytic in an extended physical sheet. Define 

F ai ,a 2 ;n(9) := F^ ai l ^^ a2 ^ n {9). (19) 

This is analytic in the extended physical sheet Im($) E [0, 27m] except for possible bound state 
poles and for kinematic poles with specific residues. It takes the values of in-state matrix 
elements on (0,oo) and {2-Kin — oo,2mn), and of out-state matrix elements on (—oo,0) and 
{2mn, 2mn + oo); and on the extended physical sheet the relation 


-^Q!2,Q!i;n(27ri?^ 9 ) — F ai ^ a2 - n {9) 


( 20 ) 


holds. The function is analytic except for possible bound state poles on Re(0) = 0, lm(0) E 
(0,7r) U {2ku — k, 27m), and except for kinematic poles at 9 = in and 9 = 2k in — iK with residues 
i and — i respectively. Further, within the extended physical sheet, on strips of widths 27tz it 
takes the values of form factors with shifted copy numbers. Specifically, the relations 

F^ lt fj, 2 -n(^) = F ait oi 2 . n {2Ki{ji ~ j 2 }ra + 9), (21) 

where p\ = (ai,ji) and p 2 = (g 7 >J 2 ), hold for all ji,j 2 £ {l,...,n}; we dehne {j} n as the 
smallest non-negative integer in the set { j +nk : k E Z} (that is, here we have { j\ —j 2 } n = 3 1 — 32 
if j 1 > 32 and {j\ — j 2 } n = j l — 32 + n if ji < 32 )- Note also the trivial symmetry under cyclic 
permutation of the copies, 

F(ai,j),(a2,k);n(@) ~ ^(ai,{j+l} n ),(c*2,{k+l} n ),n(^)- 

Relations (20) and (21), as well as the residue and boundary conditions, are consequences of the 
following equations, which hold in any massive QFT: 


F J 1 ,H 2 ;n( 9 l+ 2 Ki, 9 2 ) 
Fl^, n {9iM 
Res#, = e 2 Fp^ i/i2;n (0i + Ki, 9 2 ) 
Rese 1= g 2 FT uii2+1 . n (9i + Ki, 9 2 ) 


F^ 1+ i. n (02,9i) 
F^nuni^Oi) (jl/j'2) 
* °jl,j2°<Xl,Ol2 
~* Vjl ,72^01 ,0.2 


(23) 


where p denotes the anti-particle {a, j) and p + 1 denotes the copy-shifted particle ( a,j + 1). 

Additional poles are present whenever there are bound states. We will not need many details 
of these singularities, except for the following. Whenever a bound state of the particles ot\ and 
a 2 is present in the spectrum, with particle a\ coming from the left and a 2 from the right, 
the function F ait0l2 . n {9) possesses a pole at a purely imaginary rapidity. Such a pole is on the 
lower band of the extended physical sheet, at a position Kiy with y E (0,1), and has residue 
*7ai,a 2 ;n(y) ^ C proportional to the one-particle form factor of the branch-point twist field. 
Given aq and a 2 , the set of values of y corresponding to bound state poles on the lower band 
will be denoted B ait0l2 . By the relation (20), similar poles occur on the upper band for all 
bound states corresponding to the particles aq, a 2 being in the opposite order. These are at the 
positions 2Kin — Kiy and have residues —*7a 2 ,ai;n(y) for all y E B a2m . Besides bound state and 
kinematic poles, there are no other singularities on the extended physical sheet. Note that many 
models posses parity invariance; in this case, B ai ^ a2 = B a2m , and bound state poles occur in 
pairs with opposite residues, one on the lower band and one on the upper band. However we 
will not need to assume parity invariance. 
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For convenience we will also use the notation 7 Ql ,a 2; n(l) = <5«i,a 2 f° r the kinematic residues 
at y = 1, present both on the lower and upper band. We recall that the mass m! of the bound 
state y £ B Ql ^ a2 of particles aq and ct2 of the same mass m is given by m' = 2mcos(7ry/2). 

Note that the hermitian conjugate of the twist field is simply T t = T, associated with the 
opposite cyclic copy permutation, 

n (vac I T| 0], 0 2 ) jUi ,H 2 ;n = n{vac\T\9 1 0 2 ) (24) 

where —y is interpreted as (a,n + 1 — j). Another relation that will be helpful is crossing 
symmetry, 


<020i|T|vac>„ = n(vac|T|6h0 2 ) (25) 

Clearly, the left-hand side is the complex conjugate of a form factor of T\ and combining 
crossing symmetry with (24), we obtain 

tfunnW 1’ 0 2) = (^ lHai ;„(<Ml))* ■ (26) 

and in particular 


K lW n( e ) = F *2,on-,n(-0)- (27) 

Equation (27) implies that the set of bound state poles of Fa 2 .a\ (9) is the same as that of 

F ai , a2 (0 ); 


B 


Ot 1,02 


Bo.2,ol\ 


(28) 


and that the residues are related to each other by complex conjugation: 


7*1,02(3/) = 7 a 2 ,ai(y) V V e ,* 2 - ( 29 ) 

In particular, 7 Q]S is real for all a. 

We note that, due to relativistic invariance and invariance under permutation of copies, 
one-particle matrix elements F^ a .y n {0) are independent of 9 and of the copy number j; we will 
denote them by F a - n . By similar arguments as those above, the following relation holds: 

(■ F a -n )* = Fa-n- (30) 


Finally, in the calculation below we will need to perform an analytic continuation in the 
values of n. This analytic continuation is ambiguous, and usually, for twist field form factors in 
integrable models, one lifts the ambiguity by demanding that the large-n behaviour be described 
by a series of increasing powers of n -1 . In the present case, since form factors are generic, this is 
hard to verify. We will instead assume that matrix elements of branch-point twist fields where 
all sheet numbers are equal to 1 have analytic expressions with correct behaviour at large n. 
We will assume that this analytic expression is such that the positions of poles in the lower part 
of the extended physical strip do not depend on n, and that of poles in the upper part depend 
on n linearly, as 2nin — iriy (so that in particular the set B aim does not depend on n). We 
will also assume that these form factors vanish exponentially at large rapidities on the extended 
physical strip, and that they are not singular at n = 1/2. All these properties will be verified 
explicitly in integrable models [41]. 


4 Universal scaling of the negativity 


The results (12) and (13) indicate that a computation of the negativities and will require 
the evaluation of the two point functions n (vac|T(0)T 2 (r)|vac) n and n (vac|T(0)T(r)|vac)„ for 
n even. Such computations are routinely performed by inserting a decomposition of the identity 
in terms of asymptotic states 

1 = ^ ^ f (2 Xfc, ! I^ 1 • • • Qk)n' ' ' 0k\, (31) 

k=o (ib-A y > 

between the two fields involved. In the expression (31), the integral over rapidities is unrestricted: 
this includes the basis of in-states and out-states, as well as generalized bases obtained by 
analytic continuation to other rapidity orders. The factor of kl in the denominator accounts for 
this overcounting. Upon this insertion, the twist field correlators (and their logarithms) may 
be expressed as infinite sums over products of the form factors of the operators involved. It 
is well known that this sum naturally provides a large-distance expansion, where terms with 
higher numbers of particles decay faster as the distance r increases. Here we are interested 
in the leading (saturation) and next-to-leading order contributions to the negativity and this 
means that we will only consider contributions up to two-particle states. Let us now analyse 
each correlator in more detail. 


4.1 Two adjacent regions 

Let us take the expression (12) as our starting point. First, the behaviour at short distances 
can be obtained by using the zeroth order of conformal perturbation theory: 

ra (vac|T(0)T 2 (r)|vac) n ~ C n r _4A ™/ 2 n (vac|T(0)|vac) n (n even) (32) 

where C n is the CFT the structure constant fU; as shown in [30] (see Ecp (90)), this equals 

. Here we used the fact that the conformal dimension of T 2 is given by (11) when n is 
even. Using the definition (7), analytically continuing to n = 1 (from even values of n) and 
subtracting log(Ci) (through a re-definition of the non-universal short-distance cutoff e), this 
gives the first line of (8). On the other hand, in the limit where the distance r is infinitely large, 
the two-point function factorizes and only the zero-particle term contributes. Thus (for n even 
and using (10)) the function factorizes into vacuum expectation values: 

n (vac|T(0)T 2 (r)|vac) n -> n (vac|T|vac) n « (vac|T|vac)|. (33) 

Since i(vac|T|vac)i = 1, employing the definition (7) and again subtracting log(Ci) this implies 
that the negativity saturates to the value (9). The value of £ sat is model dependent and generally 
hard to evaluate. Expressions may be obtained for the vacuum expectation value of the twist 
fields in free models as reported in Appendix A, and we will provide a numerical analysis of the 
structure constants in a future work [41]. 

Let us now consider the next-to-leading order correction to (9). This can be obtained in 
all generality by analyzing the contribution of the two-particle terms in the spectral expansion 
of n (vac|T(ri)T 2 (?’ 2 )|vac) n . This is because the contribution from one-particle form factors 
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will always vanish when taking the limit n —> 1. Written in terms of form factors, and after 
performing one rapidity integral and expanding the logarithm function, one finds 

^ =-| 1 og(mie)+<?sat+lini — -2 ^ / dd Fj UIJta . n (0) \Fj UIJt2 . n {0)) K 0 (M auaa (6) r)+... 

n even ^ ' /i,i ,/X2 00 

(•5D 

where here and below the ellipsis represents higher-particle contributions, and where 


M, 


01,02 


( 0 ) = / 


m Ol + m Q 2 + 2m c 


, TTirvo cosh 9 


(35) 


is the square-root of Mandelstam’s s-variable. We will show that this leading correction is 
exactly 


£ 1 = - 7 log(mie) + £ S at - ~ ^ y, K 0 (V 3 m o r) 


3\/3 


ot y£l 

n(i, 


7 T 


cosec 


“ y£B at a 

n(i/2,i) 

p0 ° 7/1 /tt-z 0 
„ n / i d9 cot —-— 
8vr 2 ^ V 4 


7ry + 7T 


A'o I 2m„ cos 


27ry — 7r 


6 


^ , 7ri + 30 . _ ^ 

1 2 j _ a; 5 “il 


Ko(2m a cosh(0/2)?’) 


+ .... 

(36) 

The higher-particle contributions are difficult to analyze in all generality, but an analysis in the 
massive free boson model [41] shows that they are 0(e~ Zmir ) for Z > 2, and we expect this to 
hold in general with mi the smallest mass of the spectrum. 


Equation (36) is the main result of this section, and is valid for any massive unitary one¬ 
dimensional QFT, integrable or not. The terms on the first line are the dominant terms, and 
give the second line of (8). Note the striking fact that, although the two-particle contribution 
to the twist field correlation function is of order 0(e~ 2mir ) for any integer n, the analytic 
continuation in n changes the order of the leading correction to saturation into a slower decay 
given by the above 0(e~^ mir ). We also remark that the discrete contributions from bound 
states (the second line in (36)) only come from the least massive “half” of the possible bound 
states between a particle a and its antiparticle, as y £ B(a, a) n (1/2,1) implies masses in the 
range (0, \/2m a ). These contributions give exponential decays with exponents between \/3m a 
and 2 m a (ordered from the least massive to the most massive allowed bound state): 


V3m a < 2 m 0 


cos 


2ny — it 
6 


< 2 m c 


(37) 


The integral contribution (the third line of (36)) is a sum over a of terms which are 0(e 2m “ r ). 
We note that it is real as 


F a-a;' 2 


7 xi + 39 


F a,o:, ’ 




1 

2 


by the crossing properties (27) and (30). Further, if every anti-particle is the particle itself, then 
one can symmetrize the cotangent factor by making the replacement 


/ 7r — i,9\ 1 1 — tanh 2 (0/4) 

V 4 J 2 1 + tanh 2 (0/4) ’ 
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as all other factors are symmetric, in particular thanks to property ( 20 ). 


Result (36) is shown as follows. Recall that the next-to-leading correction for the negativity 
comes from the expression (see (34)) 


1 f°° 

dBf aiia2;n (0)K o (M auaa (d)r) 

^ > a u a2 J -°° 

where the function f ait a 2 ;n{ 9 ) is given by 

n 

/««;»(»):= £ F Ju n J»){ F n,K,r. 


(38) 


31,3 2=1 


(39) 


We will now consider n even for the rest of this section. The factorization property (10) implies 
that the double sum naturally decomposes into sums over even and odd indices: 


/«i,a 2 ;n(0) = + /a£« 2 ;n(0) + fa’^rA 6 ) + f«tarA 9 ) 


1,0 


fid 


(40) 


where fal]a 2 2 -,n( 9 ) denotes the restriction of the sum in (39) to summation variables ji,j 2 such 
that j\ = k\ mod 2 and j 2 = k 2 mod 2. Let us hrst concentrate on fa®a 2 ;n( 9 ) and fA,a 2 -,n(fi)- 
One observes that these are identical using the relations (10), ( 22 ) and ( 21 ): 


31,32 even 


£ 4L« i „(0>R T 

jl,32 odd 


Ml + l 

2 ’ 2 ’ 2 


£ F k, K ;n(°) { F Jj-,^ 

I ” 1 

^2 F ai ,a 2 -,n( A7r ij + 0) [F ai)0l2 .^(2'Kij + <9)) 


(41) 


n 

2 


3=0 


where the algebraic operations on /x = (a,j) are interpreted as operations on the copy number 
j. Adding both and using crossing symmetry (27), we obtain 

n/ 2-1 

fa°a 2 ;n( 0 ) + faiwA 9 ) = n i? « 1 ,« 2 ;n(4vrij + 0)F dl2}Bivn/2 (2mj - 0). (42) 

3=0 


In order to analytically continue the sum in n, we separate the term with j = 0. This term is 
a product of form factors whose arguments stay within the extended physical sheet for all n > 1, 
hence we may directly analytically continue the resulting integral. Clearly T Q1)Q2; i is zero as it 
is the two-particle form factor of the identity field, and no pole of F aitQl2 - n cross the integration 
contour (the real line) as the limit n —> 1 is taken. The function F^ ^-i /2 is non-zero, and 
the kinematic pole of F^ 2 ^ rn /2 h 1 the upper band is brought to the real line from above upon 
taking n —> 1 + , as 2iti(n/2 ) — in = 0 for n = 1. Under the integration in (38), the analytic 
continuation is that where the contour avoids this pole from below. The function F &2t5iy _ i / 2 
is otherwise expected to be finite (which we have explicitly checked in free models and in the 
integrable sinh-Gordon model). Hence, with T QliQ2; i = 0, the result of the integral is zero at 
n = 1, and we conclude that this term does not contribute to the negativity. 

The sum from j = 1 to j = n/2 — 1 is evaluated by a standard use of Cauchy’s theorem. Let 
s(z) := F Qlt0l2 ; n (4:iriz+9)Fa 2t a 1 ;n/2(2niz—9) 1 which is regular at the points z E { 1 , 2 ,... , n/ 2 — 1 }. 


11 





Cauchy’s theorem gives the following expression for the sum of s(j): 



where C is a counter-clockwise coutour that surrounds the points in {1, 2 ,..., n/2 — 1} but not 
those in Z \ {1, 2,..., n/2 — 1}, and where {zf, z 2 , - - -} is the set of poles of s(z) surrounded by 
C. For the present purpose, C is chosen to be composed of two vertical lines on Re(z) = 0 + and 
Re(-s) = n/2 — 0 + . The function s(z) vanishes exponentially at z —> ±ioo, so we omit the pieces 
at ±zoo that close the contour. On the vertical line Re(z) = 0 the function s(z) is a product of 
form factors with real arguments, and using (20) the same observation holds for Re(^) = n/2. 
Hence for the reasons explained above, the analytic continuation to n = 1 of these products 
vanishes, whereby for evaluating the negativity we may neglect the integral part in (43). 

Note that the integration contour stays within the extended physical sheet, so that we can 
use the known analytic properties there. Recall that B Ql)0l2 characterizes the set of bound-state 
singularities of F aiOl2 . n / 2 (0) (if any) in the lower band of the physical strip, and that B a2}0ll 
characterizes that in the upper band, and recall the relation (28). By the principles stated in 
section 3, the function s(z) has poles at the points z* with residues r* as follows: 


* = y_J_ . * = wM p 

1 4 Ani ' 1 4t r Q2 ’ m; 2 


niy 3 9 

~Y ~ Y 


z ‘ 2 = y - + ± : r; = 

2 2 2ni 2 


'Ja2,ai-,n/2{y) 


2n 


Fa u a 2 ;n (27 Tiy + 39) 


for each y £ B ait0l2 , and 

n y 6 
Ami 
9 


4 = 


* _ 7a2,ai;n(y) p 

' 3 - -F rvn rv 


* n y . 

^=2-2 + 2 Vi ■ rl = ~ 


4n - « 2 ,«i; 2 
r Yai,oi2-,n/2(y) 


r mn — 


niy 39 


2 7T 


2 2 
F ai ,a 2 -,n {2nin - 2niy + 39). 


for each y £ B a2t0ll . In addition, the function s(z) has kinematic poles if and only if a\ = a 2 , 
four poles given by the above values of z\ and r* with y = 1 and 'y a ,a-,niX) = 1. Recall that each 
bound state residue is proportional to a one-particle form factor of the branch-point twist field 
(independent of the rapidity and of the copy number), and note that 'y ai ,a 2 -.i(y) = 0 f° r evey 
y £ B ai>a2 , and that 7 ai , a 2 ;i/ 2 (y) is expected to be finite. The contribution of the residues to 
faua 2 ;n(9) + fai,a r ,n{0) sums U P to 


n 

A 


y, 7a 1 ,a 2 ;n(y ) COt f ^ + — 1 F^ ^ 


y&B aiiC , 2 u{i} 


0^2,Oil' 2 


niy — 3 6 


+ 27 , 52 ,( 5 ! ; f(y) cot ( ?ry o l0 ) F aua2 . n (2niy + 39) 


n 

A 


E l' ny — i9 

7a 2 ,oi;n(y)cot ( -;- I T7, s 


y&B a2iCll u{l} 




niy + 39 


(44) 


ny -J- iO 

+ 27a!, 02 ;^ (y) cot ( -) F a2>ai . n (2niy - 39) 
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where we used the fact that n is even in order to simplify some of the trigonometric factors. 
Note that the sum on y runs over the bound state poles y G B a 1)Q2 and y G B a2)Dll , and the 
kinematic pole y = 1 (where all 7 ’s specialize to 1 ). 

The contribution of (44) to the integral (38) is then 


n 
87r 2 


EE/ 

ai,a 2 yeB ai ,c 2 U{1} 00 


X 


cot 

+ 2 cot 


dO K o (M ai>aa (0)r) x 
7T y — id 


4 

7r y — id 
2 


r yai,a2' 1 n(.y)B a20ll -i 


iriy + 3d 


(45) 


7a 2 ,a i; f (y) F ai,a2-,n( 27ri y + 36 0 


where we have used the fact that the integrand is invariant under d —> —d, and we have use a 
change of summation indices aq, a 2 1-7 a 2 , a 1 in order to have a single y sum. The n-dependence 
is now explicit and analytic, and since the arguments of every form factor involved stay within 
their extended physical sheet for all values of n > 1 , we may directly take the limit n —> 1 in 
(45). 

Substituting n = 1 in the integrand of (45), the only nonzero terms are those corresponding 
to the kinematic residue contribution in the first line inside the square brackets, with aq = ci 2 
and y = 'y n = 1 (this is because the form factors and bound-state 7 ’s at n = 1 are zero). This 
gives 

“ 8/ dd Ko(2m a cosh(#/2)r) cot F a ,a± ' ( 46 ) 

The limit n —> 1 of the analytic function (45), taken from large enough values of n, receives 
additional contributions: the residues of the poles that cross the integration contour as the limit 
is taken. 


Consider the first line inside the square brackets in (45). Because of the factor 7 aiiQ!2;n , the 
only nonzero contributions as n —> 1 come from the kinematic pole: we restrict the particle- 
species sum to aq = d 2 =: a, and the pole sum to the single term with y = 1 (at which the j’s are 
equal to 1). In this case, the form factor in the parentheses specializes to -E Q a ;n / 2 (( 7 rz + 30)/2), 
which has argument on the line with imaginary part 7 r/ 2 . Recall that its poles (kinematic 
and bound-state) on the upper band of the extended physical sheet are at 27ri(n/2) — iriy' for 
y’ G (0, 1] (the bound-state values y’ G Ba, a and the kinematic value y 1 = 1). As n —» 1, the 
form factor goes from a high-copy-number to copy number 1/2. As n reaches 1, the poles on the 
upper band have reached the segment [0, iri). The set of these poles in [0, Tri/2 ) have crossed the 
integration line; these are for 1 / G (1/2,1]. The residues taken have 0- value, in the limit n —> 1, 
given by 


7 ri + 3d 
2 


m(l — y') =>■ d 


mil - 2y') 
3 


and the part of the residue coming from the form factor itself, F a dt . n / 2 {{Tri + 3d)/2) is, in the 
limit n -» 1, given by — 2 / 7 a, a; i/ 2 (y / )/ 3 . Evaluating the corresponding Cauchy contributions 
on the first line inside the square brackets in (45), with the factor 27ri as the Cauchy contour 
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around each pole is counter-clockwise, we find 

K 0 {V3m a r ) - ^.aji (2/) tan 

Q! J/€Ba,a 

n(i/2,i) 

(47) 

where the first term is for the kinematic pole (y = 1) and the second for the allowed bound state 
poles as described (under the change of summation index a4a). 

Now consider the second line inside the square brackets in (45). There, the form factor 
F ai ,a2-,n{2,niy + 30) has argument on the line with imaginary part 2iry. Recall that its poles on 
the upper band of the extended physical sheet are at ‘Inin — niy' for y' e (0,1] (bound-state 
and kinematic values). The subset of these poles satisfying 2n — ny 1 < 2ny cross the integration 
line upon taking the limit n —> 1. Note that for y < 1/2, none of the poles cross. The bound 
state poles have residues 7a 2 ,a i; n(j/)) which vanish at n = 1. Hence, the only pole that will give 
a contribution is the kinematic pole, at y’ = 1, which crosses the integration contour if and only 
if y > 1/2. Hence, again this constrains the particle-species sum in (45) to a\ = ct 2 =: a and 
the sum over bound state / kinematic positions to y > 1/2. In the limit n —> 1, the 0-value of 
the crossed kinematic pole is 



ny + 7r 
6 


it 0 I 2 m n cos 


2ny — 7r 
6 


2niy + 39 = ni =7- 9 = 

and the part of the residue coming from the factor F aj a-n{2niy + 39) is —i/3. Evaluating the 
corresponding Cauchy contributions on the second line inside the square brackets in (45), again 
with the factor 27ri as the Cauchy contour around each pole is counter-clockwise, we find 


K 0 {V3m a r)-~ ^ 7 Q ,a ; §G/) cot 

a ySB a ^ 
n(i/2,i) 

(48) 

where the first term is for the kinematic pole y = 1 and the second for the allowed bound state 
poles. 

Finally, we analyze /afU 2 ;n(0) and fa’i,a 2 ;n{9), the sums over mixed (even-odd) indices. In 
this case, the form factors of T 2 factorize into products of one-particle form factors (which are 
independent of the rapidity and of the copy number): 



ny + 7r 
6 


Kq I 2 m n cos 


2ny — n 
6 


fO’l (ft) = fbO 
J ai,a2',n\ u J J OL2,OLi m i n 


(-«)= E 

Jl even 
J 2 odd 


n 

2 


/* ] F ait0l2 ; n (2ni(2j + 1 ) + 0)F av nF a2 .n 

j =0 


(49) 


where we used (30) and the second equation of (23). We again use formula (43) in order to 
evaluate the sum, with s(z) = F aija2 - n (2ni(2z + 1) +9) and choosing a counterclockwise contour 
formed by the two vertical lines Re(z) = —1/2 and Re(^) = n/2 — 1/2 (these are the lines that 
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bound the extended physical strip). The result is 



where the sum is over the bound states y E B ai>ol2 (for the lower band of the extended physical 
strip) and y E -B Q2)Ql (for the upper band), and the kinematic pole y = 1 (which contributes 
if and only if a± = d^)- Putting this inside the integral in (38) and taking the limit n —> 1, 
only the term in the sum corresponding to the kinematic pole contributes, as all other terms 
gives '-fa 1 ,a 2 -,i = 0. The result, including both contributions from integrating fai,a 2 ;n(S) and 
fai,a 2 ;n{0 ) (which are the same), can be written as 

COt {^T^) K 0 ( 2m « COsh (^) r ) ■ ( 51 ) 

Combining (46), (47), (48) and (51), we find (36). 

4.2 Two semi-infinite non-adjacent regions 

Consider now the function (13). The short distance behaviour is obtained by zeroth-order 
conformal perturbation theory, 

n (vac|T(0)T(r)|vac) n ~ C n r _4An+4A "/ 2 n (vac|T 2 |vac) ri (52) 

where we used (11) for n even, and conformal normalizations of the fields. Again C n is the CFT 
structure constant Cy^-. This gives rise to the first line of (14). 

On the other hand, the leading behaviour of the logarithmic negativity at large distances is 
dictated by the factorization 

n (vac|T(0)T(r) |vac) n ->• n (vac|T|vac)^. (53) 

Thus, the negativity will decrease as r —> oo to a constant value which in this case is simply 
2 log(i(vac|T|vac)i) = 0 for all theories, in agreement with physical intuition (there should be 
no entanglement between infinitely separated regions). The entanglement build-up as regions 
get closer is obtained again by means of a form factor expansion which takes the following form 

1 f°° 

£^= lim / ^4^)0) ^(J)r) + .„, (54) 

n even fl\ ,/i2 

where ellipsis are higher-particle contributions, and where M aitOl2 {0) is as in (35); the one- 
particle form factor contribution is once again not written explicitly as it is vanishing for n —> 1 . 
We will now show that this leading correction is exactly 

£ Hh = 7^2 \ m a rK o{m a r) 2 + K 0 (m a r)Ki(m a r) - m a rKi(m a r ) 2 ] + ... (55) 
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giving the second line of (14), where the higher-particle contributions are expected to be of the 
order of the 4-particle contributions to the two-point function. 


Let us consider the sum 

n 

9a u ar,n(0) = E rfuiW (*) K, W ;n(*))’* ■ (56) 

jl J2=l 

In this case we are summing over all indices (even or odd) so this is a much easier sum to perform 
than (39). We can simply exploit the properties of form factors (21), (22) and (27), as well as 
(24), in order to write 

n— 1 

9ai,a2',n{9) = Tl ^ ( Foti,ct2',n{2TTij ~\- 9) 7^5 1) Q 2 ;n(27rij T 9). (57) 

j =0 

The sum above may be computed in much the same way as the sum (42), namely by exploit¬ 
ing Cauchy’s theorem exactly as in (43) with n/2 replaced by n and s(z) = F ai ^ a2 - n {2-Kiz + 
9)Fa 1 ,a2\n{^iz + 9). The integration contour C is again taken to be composed of two vertical 
lines on the boundaries of the extended physical sheet. As before, the part of the integration on 
the vertical lines vanish as n —> 1 , hence can be neglected for the evaluation of the negativity. 
Therefore, we concentrate on computing the residues of the poles falling within the integration 
contour. In general we will again have both kinematic and bound state poles. The kinematic 
poles are Z\ = \ — ^ and z 2 = n — \ — ^ and they occur when a\ = a 2 . The bound state 
poles are | ^ anc ^ n ~ % ~ 271 an< ^ tlmy occur for all y E B ai:OC2 U B a2jai . The bound state 
poles will necessarily give contributions proportional to the residues j ai ,a 2 ;n(y) or 7 a 2 ,ay,n(y) (or 
their complex conjugates). As these vanish as n —>• 1, we may neglect their contribution for the 
evaluation of the negativity. Hence we restrict to a\ = 07 =: a and only consider the kinematic 
poles. 


The kinematic poles of the form factors give rise to double poles in the function Fa^-ni^ij + 
9)Fa ja;n (2irij + 9 ), so the computation of the corresponding residues is a bit more involved. Let 
us consider the form factors involved. Near the kinematic poles, they have the form 


Fa,oc,n(9) 


— VK 


T F a a - n {9^ 


—1 


9 — 2irin + in 


+ Kc 


(58) 


for some constants n a , where we used the property (20), F Q!S;n (27rm — 9) = Fa t a- n (9), in order 
to relate the constants on both poles. Let us compute the corresponding residues of the function 
s(z) cot nz. Around the pole z\ we have the leading divergency 


2ni(z — z\) 


+ K, c 


2ni(z — z\) 


{na + Ka) cot nz\ cosec 2 7 rzi 
+ Ka cot nz = ... ~\ -- —7 -r--—;-r + . . . 


2n(z — z\) An(z — z\) 


and 2ni times the residue then gives 

i(tz a + Ka) cot nz\ — 


2 cosh 2 (9/2) 


(59) 

(60) 


Similarly, around Z 2 , we find 
— i 

+ Ka 


2ni(z — Z 2 ) 


— 1 \ (Ka + Ka) COt nZ 2 COSeC 2 7 TZ 2 

+ Ka cot 7 TZ = . . . --—7-r--—7-r + . . . 


2ni(z — z 2 ) 


2n(z - z 2 ) 4n(z - z 2 ) 


(61) 
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and 2iri times the residue then gives 


— i(K a + Ka) COt TTZ2 


i 

2 cosh 2 (0/2) 


(62) 


Adding these two contributions, using cot 7 rzi = cot 7 r ^2 = — tan(0/2i) and finally dividing by 
—2 i and multiply by n we get the exact contribution to the sum, 

|sech 2 ^, (63) 

from the kinematic poles. Inserting the contribution (63) into the integral (54) and taking the 
limit n->lwe obtain the following leading correction to the negativity: 

-I _ roa n / f)\ 

£ Hh = -^2 J d0 seek 2 -I< 0 (2m a cosh r). (64) 

This integral may be computed exactly giving the universal result (55). 


5 Conclusions and outlook 


In this paper we have considered the logarithmic negativity, a measure of bipartite entanglement, 
in a general 1+1-dimensional massive quantum field theory, not necessarily integrable. We have 
considered two particularly interesting limiting cases. Firstly, we have computed the negativity 
between a finite region of length r and an adjacent semi-infinite region. Secondly we have 
computed the negativity between two semi-infinite regions separated by a distance r. The main 
results of our analysis can be summarized as follows: 

1) Both limits of the negativity tend to a constant for r —> oo. In the latter case, this constant 
is zero, which agrees with physical intuition according to which there is no entanglement 
between infinitely separated regions. For adjacent regions, the saturation constant is 
theory-dependent, requiring the knowledge of the expectation values of twist fields and 
of the twist field three point coupling Cfj-j-. We have provided explicit formulae for the 
contribution to saturation of these expectation values in free Majorana and free Boson 
theories (see Appendix A). However, general formulae for the three point coupling involved 
are not yet known, even for free models. 

2) The leading corrections to saturation are (in both limits) exponential decays in r (described 
by modified Bessel functions) that are solely controlled by the mass spectrum of the model, 
independently of its scattering matrix. This implies that, like the entanglement entropy, 
the logarithmic negativity displays a very high level of universality, allowing one to extract 
information about the mass spectrum. 

3) By studying sub-leading terms we have shown that, unlike the entanglement entropy [23, 
24], a large-r analysis of the negativity of adjacent regions allows for the detection of bound 
states. 

There are still many additional features of the negativity we would like to explore for 1+1 
dimensional QFT. From the numerical point of view, it would be very interesting to test our 
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present predictions numerically by evaluating corrections to the negativity of gapped quantum 
spin chains. It should also be possible to exploit the form factor approach in order to obtain 
numerical estimates of the CFT three point coupling C^fj- in free theories. From the point of 
view of analytical computations, it would be natural and interesting to analyze further sub¬ 
leading corrections to the limiting cases of the negativity considered here. Preliminary results 

[41] suggest that some of these corrections may still be universal. For the case of two semi¬ 
infinite regions separated by a distance r we expect that a full description of all corrections may 
be accessible for free theories by performing an analysis reminiscent of the study carried out in 

[ 42 ] - 

A more general study of the negativity will involve the evaluation of more general matrix 
elements of the form n {9 \, • • • , 9k\T\j3 p ■ ■ ■ /3i) n with all the subtleties that are well known for 
this kind of problem and the additional complication of considering twist fields. We expect to 
report progress in this direction for integrable models in the near future [41], 


A Formulae for the free Majorana and free Boson theories 


The general results (36) and (55) can be easily specialized to the free Majorana (FM) and free 
Boson (FB) theories (with mass m ). In the case of (55) this is immediate whereas (36), as we will 
see, exhibits some model dependence in the subleading corrections. The FM and FB theories 
are very simple, consisting of a single particle spectrum with no bound states and two-particle 
scattering matrices given by S(6) = —1 and S(9) = 1, respectively. Due to internal symmetries, 
branch point twist field form factors are only non-vanishing for even particle numbers (even for 
general n ). In particular, the two-particle form factors associated to two particles on the same 
copy (normalized by the vacuum expectation value) are given by: 


-iFM 


(*) = 


sin - 

n 


sinh ± 


and F™(9) = 


2nsinh(^)sinh(^) sinh £ 


sin - 

n 


2ns inh(^)sinh(^)’ 


(65) 


and all higher particle form factors may be obtained by employing Wick’s theorem. For the 
FM theory these form factors can be written in terms of a Pfaffian determinant [42]. The same 
formula may be used for the FB theory by simply replacing all terms with negative sign by 
terms with a positive sign. 


The values of £ sa t and £ s hift given by (9) and (15) are model dependent and generally hard 
to evaluate. Although we do not have expressions for the CFT structure constant involved, 
expressions may be obtained for the twist field vacuum expectation values in free models, thus 
giving ^-shift+gsat (16). For the FM theory the expectation values of twist fields are known [23] 
for generic values of n and the following value can easily be obtained: 


cFM 

c shift 


I cFM 
> 


■'sat 


i 

8 


log 2 + 



dt 
2 1 


r e -2t 


+ 


1 — cosh t 
sinh t sinh(2f) 


1 log 2 3 log A 

- 4-— H-— 

8 3 2 




-0.06319... 


( 66 ) 


where A is Glaisher’s constant A = 1.28242.... In the FB theory, the branch-point twist held 
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has vacuum expectation value 


n (vac|T|vac)n 


. n-i (Tin) 
(27T)— ' 1 ' 


n“ 



n coth | — coth 


2(1 - e 





( 67 ) 


where in A„ the central charge is c = 1, and this gives 

cFB , cFB i 4 

shift sat = - log - = 0.06039... ( 68 ) 

2 4 7T 

In both cases the expectation values may been obtained by employing the angular quantization 
approach developed in [43, 44], The free Boson computation will be presented in detail elsewhere 

[41]- 

Concerning the other terms in (36) the integral correction in the third line is in fact zero 
for the FB theory as all form factors vanish in the limit n —> On the other hand, in the FM 
theory with mass m, this correction is exactly given by 


1 



dO sech-sech 
2 


36 

T 


Ko(2mr cosh(0/2)). 


(69) 


We thus have overall 


£ 1,fm = —-log(me) — 0.06319 + log(Ci) 
8 


3v/3vr 


1 9 3 q 

A'o(v / 3m.r) + -—/ dQ sech-sech— I\ 0 (2mr cosh(0/2)) + • • • (70) 

( 27r ) 4-oo 2 2 


and 


\ log(me-) + 0.06039+log(C'i)- ^^Ko(V3mr) + 

4: 3v37t 


(71) 
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